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Abstract 



We present the results of an exploratory numerical study of two dimensional QCD with 
overlap fermions. We have performed extensive simulations for U{Nc) and SU (Nc) 
color groups with Nc = 2,3,4 and coupling constants chosen to satisfy the 't Hooft 
condition g'^Nc = const = 4/3. We have computed the meson spectrum and decay 
constants, the topological susceptibility and the chiral condensate. For U{Nc) gauge 
groups, our results indicate that the Witten-Veneziano relation is satisfied within our 
statistical errors and that the chiral condensate for Nf = 1 is compatible with a non- 
zero value. Our results exhibit universality in Nc and confirm once more the excellent 
chiral properties of the overlap-Dirac operator. 



1 Introduction 



Several years ago, in a pioneering investigation, 't Hooft studied U{Nc) gauge theories in 
the limit Nc —>■ oo with Ct = g'^N^ kept constant He showed that only planar diagrams 
with quarks at the edges dominate and therefore some non-perturbative QCD physical 
observables can be computed in this limit. He proposed two dimensional models with 
important features of QCD, but simple enough to sum explicitly all planar diagrams in the 
meson spectrum computation. 

Recently, the overlap formulation |^ ^ has made it possible to introduce a lattice Dirac 
operator D which preserves a lattice form of chiral symmetry at finite cut-off As a 
consequence the Ua{^) chiral anomaly is recovered a la Fujikawa and the Dirac Operator 
has exact chiral zero modes for topologically non-trivial background configurations 0. A 
precise and unambiguous implementation of the Witten-Veneziano formula can be obtained 
i- 

These theoretical developments generated a renewed interest in 't Hooft's results and 
prompted us to perform an exploratory numerical investigation of a class of two dimensional 
non-Abelian models with overlap lattice fermions. Precisely, we have simulated models of 
QCD2 with U{Nc) and SU{N^) color groups for A^'^ = 2,3,4 {Nf = 0,1,2) imposing the 
't Hooft's condition Ct = g'^Nc = const. Our systems are small enough that we could 
compute the fermionic propagator and det{D) exactly following the scheme used in 



(see Refs. |P, |10| and ||TT], [T^, |13|, [Tj] for more refined implementations in two and four 
dimensions) . 

For U {Nc) models we have found many background gauge configurations with zero modes 
in the fermionic operator. By counting them and averaging over the configurations we have 
computed the quenched topological susceptibility obtaining values in very good agreement 
with the analytic results. We have computed the chiral condensate for Nf = 1 which turns 
out to be compatible with a non-zero value. For SU{Nc) models we have not found any 
configuration with exact zero modes as expected since these models have an exact Ua{1) 
symmetry in the chiral limit. 

From two-point correlation functions of fermion bilinears we have extracted the meson 
masses and the corresponding decay constants. In the U{Nc) case the t]' mass in the chiral 
limit verifies the Witten-Veneziano relation |jT5|, [1^ within errors for each Nc- The pion 



masses verify quite well the expected functional dependence oc M. For SU{Nc 



models our data favor the functional dependence oc m^' [p!7| , p!q] . In both cases 



at fixed Nj, data exhibit universality in Nc and quenched results get closer and closer to 
unquenched ones when N^. increases. 

In the next section we shall briefly remind the reader of some properties of the 't Hooft 
model in the continuum. In section H we define the overlap regularization we have imple- 
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merited numerically. In section ^ we present our numerical results for the meson spectra and 
decay constants. In Section |^ we compute the topological susceptibility and we compare the 
r/' mass with the one extracted from the Witten-Veneziano formula. In section ^ we report 
our results for the chiral condensates. Section |^ is devoted to some concluding remarks. 



2 The 't Hooft Model in the Continuum 

We consider two-dimensional models with color group U {Nc) and SU {Nc) with Nf degenerate 
flavors defined by the action 



S = d'x 



(1) 



The if) is a, A^j-dimensional fermion multipletQ and we use the following representation of the 
two-dimensional 7-matrices 

7i = (Ti , 72 = CX2 , 75 = -«7i72 = 0-3 (2) 

where cxj are the Pauh matrices. In Eq. (Q) = + igA^^, where = A^t^ is the gauge 
potential and are the iV? or N^ — l generators for the groups U{Nc) or SU{Nc) respectively, 
normalized according to ti (t^t^) = 1/2 6^^. The field strength reads F^^, = F^^t^^ where 

F^. = - d^A^^ - gf^^^'A^^A^^ (3) 

The models described in Eq. are super-renormalizable and therefore g and rriq are finite 
bare parameters, 't Hooft studied the U{Nc) models in the limit 

iVc ^ 00 , g'^Nc = Ct = constant (4) 

which corresponds to take only planar diagrams with no fermion loops P[. 



2.1 The U{Nc) models 

The massless action in Eq. (|1]) has an Uv{Nf) UA{,Nf) flavor symmetry. The UyiX) 
symmetry is preserved while the Ua{^) is softly broken by the quark mass and explicitly 
broken by the anomaly which in two dimensions appears in two-point functions, not in 
triangle loops [0. The corresponding singlet Ward identities are 

d^y^{x) = (5) 
d^A^{x) = 2NfQ{x)+2m,P{x) (6) 
^ Color and spinor indices are suppressed throughout the paper. 
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where 

Vf,{x) = ip{x)-ff,ip{x) , Af,{x) = -ie^yVy{x) = V'(x)7^75V'(x) , P{x) = 'ip{x)-f5ij{x) (7) 
and the topological charge density reads 



Qi^) = ^ e^uF'^ui^) (8) 



where is the Abelian field strength. The Ward identities associated to the non-singlet 
(A'^^ > 1) axial and vector symmetries are given by 

d.Vfix) = (9) 
d^Alix) = 2m,pf{x) (10) 

where V^, and are defined as in Eq.(|^), but with the insertion of the proper fiavor 
generator. 

The topological charge of a U{Nc) background gauge configuration is 

Q = Jd'x Qix) = J d'x t,,Fl^{x) (11) 

and it is related to the difference between the number of positive (n+) and negative {nJ) 
eigenvalues of the Dirac operator through the Atiyah-Singer theorem 

Q = n_-n+ (12) 



In two dimensions the vanishing theorem ensures that only or n_ is non-zero [T^. The 
topological susceptibility x i^i the pure gauge theory is EC 



X = {Q{x)Q{mYM = ^, (13) 

where the expectation values (. . .)\ym in Eq- (0) have been taken in the theory without 
fermion fields. 

Since the Ua{^) symmetry is anomalous, for Nj = 1 one can have a non-zero chiral 



condensate 0, ^ (V^^) without violating the Mermin- Wagner theorem It is related to 
the topological susceptibility in the full theory as 

X = -mg(V^^) + 0(mJ) (14) 

On the contrary, for Nf > 1 a non-zero condensate would spontaneously break the SUA{Nf) 
symmetry. 



The meson spectrum of the two-flavor U{Nc) models exhibits a pseudoscalar flavor-singlet 
excitation [t]') and flavor-triplet quark- ant iquark bound states (pions). The rj' is massive due 
to the anomaly and its mass in the chiral limit is 

' 71 I\r 71 



The Witten-Veneziano formula [|T^, for the U{Nc) models gives 

= (16) 

where /,r is the pion decay constant. By inverting Eq.(|T6|), assuming that it is an exact 
relation between the r]' mass and the topological susceptibility and taking into account 
Eqs. (0), (0) and (|TB[) one gets 



In the Nc ^ limit and for rUq <^ yCt/vr, the pion mass M^^ as a function of the quark 
mass rrig reads 0] 

3 

It is interesting to note that is linear at the leading order in the quark mass like in four- 
dimensional QCD and the coefficient in front of rriq is expected to be exact in the Nc oo 
limit 01. 



2.2 The SU{Nc) models 

Analogously to the previous case, the massless action (|1|) of two dimensional SU{Nc) models 
has an exact Uv{Nf) UA{Nf) flavor symmetry. Since in this case there is no Abelian field 
strength component F^^,, the Ua{^) symmetry is only softly broken by the quark mass rriq. 
Therefore the Dirac operator should not have any zero-modes. 



In the limit Nc ^ oo, a non-zero chiral condensate was obtained in 
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and a Berezinski-Kosterlitz-Thouless phase transition was advocated to reconcile this result 
with the Mermin- Wagner theorem |22|. This behavior would favor the argument that in the 
limit Nc ^ oo both U{Nc) and SU{Nc) gauge groups should describe the same physics [|T], 
23, E5i. 
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The meson spectrum of the SU{Nc) models exhibits only pions. For Nf = 1 and the 



SU{2) color group their mass has been computed in |]T8l in the semiclassical WKB approx- 
imation 
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1 ,e\4 



n 77 



(20) 



3 The Overlap Dirac operator 



We have implemented the lattice action 



S = SGiU)+ Y.^^{x)D^^{x,y)^,{y) 

i=l,Nf x,y 



(21) 



where Sg{U) is the standard Wilson gauge action 

1 



Sciu) =/3y: 



2Nr 



(3 = 2Nc/{gaY, a and g being the lattice spacing and bare coupling constant and 



U^u{x) = U^{x)Uu{x + a(x)Ul{x + aO)Ul{x) 



In Eq. m 



DmAx,y) 



'I 



TUqa, 



D{x,y) + ruq S: 



xy 



where the Neuberger-Dirac operator is defined as 



D = -{l + {Dw ^ 



a a 



Dw is the Wilson-Dirac operator 

Dw ■ 

where 

V^^i{x) = 
V>.(x) = 



7M(v^ + v;)--av;^v^ 



- [Uf,{x)ijji{x + afi) - i)i{x)] 

— ipi{x) — UUx — a(i)ipi{x — afi) 



(22) 

(23) 
(24) 

(25) 

(26) 

(27) 
(28) 



The Neuberger Dirac operator is 75-hermitian, i.e. Z)^ = 75D75, and satisfies the Ginsparg- 
Wilson relation 

75/^-^ + Z^-Ss = a75 (29) 



which guarantees that in the chiral hmit the lattice action (|2T|) is invariant under the con- 
tinuum symmetry P| 

5tpi = 75(1 - aD)'ipi, 5'ipi = ipa^ (30) 

As a consequence the Ua{^) anomaly, if present, is recovered a la Fujikawa and the 
Dirac operator has exact chiral zero modes for topologically non-trivial gauge field configu- 
rations [^. The analogous flavor non-singlet chiral transformations are obtained by including 
a flavor group generator in Eq. (^0|). 

From Eq. (^) one can derive the following identities (a = 1) 




which turns out to be useful in the numerical implementation. 



4 The Pion Masses and Decay Constants 

In order to study the response of the overlap Dirac operator in the presence and absence of 
the chiral anomaly and at the same time to analyze the scaling of the physical observables 
with Nc (at fixed A^/), we have performed extensive simulations of U{Nc) and SU{Nc) models 
for Nc = 2,3,4. We have generated the gauge configurations with a standard Metropolis 
Monte Carlo algorithm according to the gauge action in Eq. (0). To avoid the Gross- 
Witten phase transition |^ we have chosen Ct = g'^N^ = 4/3, which corresponds to (3 = 6 
for Nc = 2, (3 = 13.5 for Nc = 3 and /3 = 24 for Nc = 4. We have generated 500 independent 
configurations for Nc = 2, 300 for Nc = 3 and 150 for Nc = 4: separated by 10000 sweeps of 
the whole lattice. For all gauge groups we have fixed the same dimensions in lattice units, 
i.e. iVt = A^^ = 18. Consequently the Neuberger-Dirac operator is a complex matrix of 
dimension 1296 x 1296 for Nc = 2, 1944 x 1944 for A'^ = 3 and 2592 x 2592 for A'^ = 4. We 
could diagonalize exactly the hermitian operator D + D'^ hj using full matrix algebra routines 
with the resources available to us. By using Eq. (|31|) we have computed the eigenvalues of 
the massive Neuberger operator and its determinant and with Eq. (^) we have determined 
the propagators of the massive fermions. We have explicitly checked that, for each of the 
fermionic masses rrig = 0.04,0.05,0.06,0.07,0.08,0.1, the lattice spans at least four pion 
correlation lengths. The effects of dynamical fermions have been included by weighting the 
observables with the appropriate powers of the fermion determinant. The smallness of the 
18 x 18 lattice warrants this procedure |^ , which would lead to an unacceptable large variance 
on larger systems. 
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The meson masses and the decay constants have been extracted in the standard manner 
from the vector correlator at zero momentum 

Y,(y,{x,y)V,{x\y + t)) (33) 

x,x' ,y 

where the non-singlet vector current we have used is 

V^{x,y) = ^i(x,i/)7^[(l - '^D)^2]{x,y) (34) 
In Table |l] we report the pion mass squared for the U{Nc) models with two flavors of 





iVe = 2 


iVe = 3 


iVc = 4 


mq/gy/W^ 






Ml/g'N, 




Nf 


= 2 




0.0346 


0.046(10) 


0.047(3) 


0.055(6) 


0.0433 


0.056(10) 


0.059(3) 


0.067(5) 


0.0519 


0.068(8) 


0.072(3) 


0.080(4) 


0.0606 


0.080(7) 


0.085(3) 


0.093(3) 


0.0692 


0.094(6) 


0.098(3) 


0.107(3) 


0.0866 


0.124(4) 


0.128(3) 


0.138(2) 


quenched {Nf = 0) 


0.0346 


0.069(1) 


0.063(1) 


0.057(1) 


0.0433 


0.083(1) 


0.077(1) 


0.071(1) 


0.0519 


0.098(1) 


0.091(1) 


0.086(1) 


0.0606 


0.113(1) 


0.106(1) 


0.100(1) 


0.0692 


0.129(1) 


0.121(1) 


0.116(1) 


0.0866 


0.161(1) 


0.153(1) 


0.148(1) 



Table 1: M^/g^Nc vs. mjg^^, U{N^) models. 



dynamical fermions and in the quenched approximation. The quenched results get closer to 
the unquenched two-flavor ones when Nc gets larger. Figure ^ shows as a function of 
the quark mass for Nf = 2 and provides evidence of universality in Nc for the pion masses. 
According to Eq. (0), should have a linear dependence in and it is expected to 
vanish in the chiral limit. From a fit 

g^Nc gVNc 
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0.15 



0.125 



0.1 



0.075 



0.05 



0.025 



o Nc=2 
□ Nc=3 
A Nc=4 



0.02 



0.04 



0.06 



0.08 



0.1 



m7g(N,; 



Figure 1: M^/g'^Nc vs. rriq/ g^/Nl for Nf = 2, U{Nc) models. 



we obtained A = -0.013(17) and B = 1.57(19) for Nc = 2, A = -0.007(4) and B = 1.54(4) 
for Nc = 3 and A = —0.006(8) and B = 1.66(9) for Nc = 4. In order to carefully determine 
the systematic error affecting B, further simulations would be necessary that go beyond the 
scope of the present investigation. 

The pion decay constants /tt turns out to be almost constant in niq within our statistical 
errors. Their values in the chiral limit are reported in Table 0. 



Nc = 2 

u 


iV, = 3 N, = 4 

fn fn 


Nf = 2 


0.7(1) 


0.94(5) 1.10(1) 


quenched {Nf = 0) 


0.77(1) 


0.96(1) 1.12(1) 



Table 2: U{Nc) models. 



In Table ^ we report the pion masses for the SU (Nc) models with one and two flavor of 
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dynamical fermions and in the quenched approximation. Figure ^ shows the pion masses for 





iVc = 2 


Nc = 3 


A^c = 4 








M^/g^^ 




Nf = 


2 




0.1063 


0.183(15) 


0.207(23) 


0.181(28) 


0.1233 


0.211(12) 


0.237(23) 


0.212(26) 


0.1392 


0.236(10) 


0.263(22) 


0.240(24) 


0.1543 


0.259(8) 


0.286(19) 


0.267(23) 


0.1687 


0.280(7) 


0.308(16) 


0.292(21) 


0.1957 


0.321(5) 


0.348(10) 


0.341(18) 




Nf = 


1 




0.1063 


0.182(7) 


0.196(4) 


0.199(12) 


0.1233 


0.211(6) 


0.226(4) 


0.229(10) 


0.1392 


0.237(5) 


0.254(4) 


0.258(9) 


0.1543 


0.262(4) 


0.279(3) 


0.285(8) 


0.1687 


0.285(4) 


0.303(3) 


0.310(7) 


0.1957 


0.329(3) 


0.349(2) 


0.357(6) 


quenched (iVj = 0) 


0.1063 


0.201(3) 


0.212(1) 


0.216(1) 


0.1233 


0.227(2) 


0.241(1) 


0.246(1) 


0.1392 


0.253(2) 


0.268(1) 


0.273(1) 


0.1543 


0.277(2) 


0.293(1) 


0.299(1) 


0.1687 


0.300(2) 


0.316(1) 


0.323(1) 


0.1957 


0.344(2) 


0.361(1) 


0.368(1) 



Table 3: M^/g^/K vs. {mjg^^f'^ SU{Nc) models. 



Nf = 2 and reveals also in this case a universality in Nc already for Nc = 2, 3, 4. According 
to Eq. (pOf ), should have an m^^^ quark mass dependence and should vanish for nig = 0. 
From a fit 

A + B{—^)^ (36) 



gy/K gVK 

we got A = 0.023(26) and B = 1.52(12) for N^ = 2, A = 0.046(43) and B = 1.55(17) for 
Nc = 3 and A = —0.01(4) and B = 1.79(12) for N^ = 4, where again the errors are statistical 
only. We performed also a fit of the form 



Mjg^c = c{m/g^)j (37) 
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0.4 
0.35 

0.3 
0.25 

0.2 
0.15 

0.1 
0.05 




ONc=2 
nNc=3 
A Nc=4 



0.04 0.08 0.12 0.16 0.2 0.24 

[myg{NJ ] 

Figure 2: M^/g^c vs. {mj gy/W^)^/^ for Nf = 2, SU{Nc) models. 



and we obtained a value of 7 compatible with 2/3 for A^^c = 2, 3, 4. 

The values of the pion decay constants /t^ in the chiral limit for the SU{Nc) models are 
reported in Table 



iVc = 2 


Nc = 3 


iVc = 4 


u 






Nf = 2 


0.64(6) 


0.92(3) 


1.21(12) 


Nf = l 


0.70(2) 


0.98(2) 


1.14(3) 


quenched (A^^ 


= 0) 


0.76(1) 


0.98(1) 


1.12(1) 



Table 4: SU{N,) models. 
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5 The Witten-Veneziano Relation and the rj Mass 



The topological charge of a given background configuration can be computed by counting the 
number of zero modes of the overlap Dirac operator. Since we have diagonalized exactly the 
overlap operator, we could compute the quenched topological susceptibility x by averaging 
the square of the topological charge, normalized with the volume, over the configurations. 
In Table ^ we report the results we have obtained for the various Nc- 



N, = 2 
X 


iVe = 3 

X 


X 


0.0258(16) 


0.0298(23) 


0.0319(37) 



Table 5: Topological susceptibility, U{Nc) models. 



Even if with increasing Nc the topological susceptibility get closer and closer to the 
analytical value x = 0.0337 given in Eq. (p!3|), more accurate studies of discretization and 
finite size effects, that goes beyond the scope of our exploratory investigation, are required 
to reach a final conclusion. In fact if the magnitude of these systematic uncertainties is 
different for the gauge groups, they would spoil any extrapolation to Nc oo. For example, 
building on a trade-off between spatial and internal degrees of freedom that prompted the 



introduction of the Nc oo single plaquette model [^], one could heuristically argue that 
the theory should be less affected by finite size effects for increasing A^^^. To try to estimate 
finite size effects in the worst case, we computed the value of the topological susceptibility of 
the U{2) model for the volume = 24 x 24 and we obtained x = 0.0316(24). The difference 
of this value with the analogous one in Table ^ gives a rough estimate of the systematic error 
induced by finite volume effects on x- 

For SU (Nc) model we directly checked that there are no zero modes for each gauge field 
configuration and therefore the topological charge is always zero. 

From the values of x reported in Table |^ and the pion decay constants at Nf = given 
in Table ^ we have computed M^, for Nf = 1 given in Table ^ (statistical errors only). It 
turns out that the WV relation is satisfied for all U{Nc) gauge groups within our overall 
uncertainties. 

The singlet pseudoscalar correlation functions are given by differences between connected 
and disconnected contributions and therefore they are noisier with respect to the non-singlet 
case. To have an independent determination of the 7]' mass for Nf = 1, we resorted the 
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Nc = 3 


A^c = 4 








Analytic 


0.159 


0.106 


0.079 


Numerical from Eg. (|16|) 


0.129(6) 


0.097(7) 


0.076(6) 


Numerical from Eq.(pO|) 


0.159(8) 


0.111(7) 


0.084(8) 



Table 6: M^, from the Witten-Veneziano relation in Eq. ([iBD and from Eq . (|40|) . 



method proposed in p8|. It exploits the quenched two-loop disconnected Eg and one- 
loop connected t\ ^*^'^P contributions to the rj' propagator 

r,(t)=/cix(P(x,t)P(0,0))L_^^^ (38) 

where the singlet pseudoscalar density is defined as 

P(x,t) = V'(a;,t)75[(l - 'i^D)i;]{x,t) (39) 



The formula we have used reads ^ ^ 

^2-loop 



Ml = 2M, lim ^ ' (40) 

|t|rJ-^°°P(t) 

The results we have obtained, reported in Tab. ^, compares remarkably well with the analytic 
result in the chiral limit. 



6 Chiral condensate 

The chiral condensate can be expressed in terms of the eigenvalues Aj = 1 -|- e'^' of the 
Neuberger-Dirac operator as (a = 1) 

1 1 ((detP r^F(Q^ cos^.)). ^^^^ 



Nfg^/K'^' 2 g^,V ((detD^J^/)^; 
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where 



F{Q,mq, cos 6i 



\Q\ 



and using Eq. 



+ 



rrir, 



E 



L — cos 9i 
+ (l-^)cos^, 



(42) 



(detD^ 



Nf\ 



n 

COS 6ij^- 



1 L 



^1 2(l + cos^0 



(43) 



For all the groups we have computed the chiral condensate by using the direct definition 
in Eq. (|4l|) . As a representative example, in Fig. |^ we compare the values of the U(2) and 
SU{2) chiral condensates for Nj = 1 computed up to mq/ gy/Wl = 0.0346 in order to limit 
the finite size effects. A chiral extrapolation of these plots hints to a finite and zero value of 
the f/(2) and SU{2) chiral condensates, respectively. 

In Fig. ^ we show the chiral condensate for all U{Nc) models with Nf = 1. They have 
been computed by using the definition in Eq. (|4l|) and the one from the Axial Ward Identity 
in Eq. (p!^ obtained by neglecting the O(m^). The comparison is interesting because the 
mass dependence of the two quantities is different. Again both the determinations point in 
the direction of non-zero condensates. An interesting improvement in the determination of 
the condensates in the chiral limit could be obtained by implementing in two dimensions a 
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m /a^'N m /aVN' 

q *- c q ^ c 

Figure 4: Nf = 1 chiral condensate. Left: from Eq. (|^). Right: from Eq. (|1^). 



finite volume technique analogous to the one proposed in four dimensions in Refs. 
Moreover the use of algorithms which generate directly dynamical configurations could help 
in reducing the fluctuations induced by the re-weighting of the observables with the fermionic 
determinant. 

7 Conclusions 

We have performed an exploratory numerical study on the lattice of two dimensional models 
defined by the gauge groups U{Nc) and SU{Nc) {Nc = 2,3,4) and Nf = 0, 1,2 degenerate 
fermions introduced by using the Neuberger-Dirac operator. Our results prove that the 
computation is feasible and it would be interesting to further pursue this line of research 
with a more detailed analysis, especially of discretization and finite size effects, that goes 
beyond the scope of our present investigation. 

We have found that within our statistical errors the pion masses verify quite well the 
expected functional dependence oc and oc m^/^ for U{Nc) and SU{Nc) models 
respectively. In both cases, at fixed Nf, data exhibit universality in Nc and quenched results 
get closer and closer to unquenched ones when Nc increases. 

As expected from an analysis of the symmetries of the models, for SU{Nc) groups we 
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have not found any background gauge configuration witfi exact zero modes in the fermionic 
operator. On the other hand, many background gauge configurations with exact zero modes 
were found for the U{Nc) models. By counting them and averaging over the configurations we 
have computed the quenched topological susceptibility obtaining values in good agreement 
with the analytic results. By using the meson decay constants extracted from the two-point 
functions, our data verify the Witten-Veneziano relation within errors for each Nc- Even if 
we could not safely extrapolate our data to the chiral limit, we have produced evidences that 
the chiral condensate for Nf = 1 is compatible with a non-zero value. 

Of course, one would be interested to know how close real QCD is to the large limit. 
In Refs. 1 3^ it was pointed out, on the basis of a study of three and four dimensional pure 



gauge theories, that even the SU{2) color group is close to SU (oo). Although our analysis 
is limited to QCD2, remarkably enough we saw that even with dynamical quarks 2, 3 and 4 
colors appear to be in the same universality class, i.e. the physical quantities are degenerate 
in Nc within the errors. 

Although our study is exploratory, the results we have obtained are very gratifying 
and indicate that it would be interesting to perform a deeper analysis in two dimensions 
and eventually extend it to four dimensions. In particular one could study within our 



computational scheme the baryons that are expected to be the QCD solitons An 
interesting improvement in the determination of the condensates in the chiral limit could 
be obtained by implementing finite volume techniques in two dimensions. Moreover the use 
of algorithms which generate directly dynamical configurations could help in reducing the 
fluctuations induced by the re-weighting of the observables with the fermionic determinant 
and therefore would allow to simulate larger lattices. 



Acknowledgments 

We warmly thank G.C Rossi for many illuminating discussions. We wish to thank A. Gonzalez- 
Arroyo and H. Neuberger for interesting discussions. This work has been supported in part 
under DOE grant DE-FG02-91ER40676. The work of F.B. is supported in part by an INFN 
Postdoctoral Fellowship. He thanks the members of the Theoretical Particle Physics Group 
at Boston University for their kind hospitality. 



References 

[1] G. 't Hooft, Nucl. Phys. B 72, 461 (1974). 



15 



[2] G. 't Hooft, Nucl. Phys. B 75 , 461 (1974); see also: G. 't Hooft, in Lectures given at 
Int. School of Subnuclear Physics, 'Ettore Majorana', Erice, Sicily, Jul 11-31, 1975. 

[3] H. Neuberger, Phys. Lett. B 417, 141 (1998); Phys. Lett. B 427, 353 (1998). 

[4] R. Narayanan and H. Neuberger, Nucl. Phys. B 443, 305 (1995); Phys. Rev. Lett. 71, 
3251 (1993); Nucl. Phys. B 412, 574 (1994); Phys. Lett. B 302, 62 (1993). 

[5] M. Liischer, Phys. Lett. B 428, 342 (1998); 

P. Hernandez, K. Jansen and M. Luscher, Nucl. Phys. B 552, 363 (1999). 

[6] K. Fujikawa, Phys. Rev. Lett. 43, 1195 (1979) and Phys. Rev. D21, 2848 (1980). 

[7] P. Hasenfratz, V. Laliena and F. Niedermayer, Phys. Lett. B 427, 125 (1998). 

[8] L. Giusti, G. C. Rossi, M. Testa and G. Veneziano, |hep-lat/0108009| . 

[9] L. Giusti, C. Hoelbling and C. Rebbi, Phys. Rev. D 64, 054501 (2001). 

[10] R. Narayanan, H. Neuberger and P. Vranas, Phys. Lett. B 353, 507 (1995); 
J. E. Kiskis and R. Narayanan, Phys. Rev. D 62, 054501 (2000); 
F. Farchioni, I. Hip and C. B. Lang, Phys. Lett. B 443, 214 (1998); 
S. Hausler and G. B. Lang, Phys. Lett. B 515, 213 (2001); 
S. Chandrasekharan, Phys. Rev. D 59, 094502 (1999); 
W. Bietenholz and I. Hip, Nucl. Phys. B 570, 423 (2000). 



[11 
[12 

[13 
[14 
[15 
[16 
[17 
[1 



R. G. Edwards, U. M. Heller and R. Narayanan, Nucl. Phys. B 540, 457 (1999); 
R. G. Edwards, U. M. Heller and R. Narayanan, Phys. Rev. D 59, 094510 (1999). 



P. Hernandez, K. Jansen and M. Luscher, [hep-lat/0007015 



P. Hernandez, K. Jansen and L. Lellouch, [hep-lat/0001008 



L. Giusti, G. Hoelbling and G. Rebbi, Phys. Rev. D 64, 114508 (2001). 
S. J. Dong, T. Draper, I. Horvath, F. X. Lee, K. F. Liu and J. B. Zhang, |hep-lat/0 108020 
E. Witten, Nucl. Phys. B 156, 269 (1979). 

G. Veneziano, Nucl. Phys. B 159, 213 (1979); Phys. Lett. B 95, 90 (1980). 
P. J. Steinhardt, Nucl. Phys. B 176, 100 (1980). 
G. J. Hamer, Nucl. Phys. B 195, 503 (1982). 



16 



[19] J. Kiskis, Phys. Rev. D 15, 2329 (1977); 

N. K. Nielsen and B. Schroer, Nucl. Phys. B 127, 493 (1977); 
M. M. Ansourian, Phys. Lett. B 70, 301 (1977). 

[20] D. J. Gross, I. R. Klebanov, A. V. Matytsin and A. V. Smilga, NucL Phys. B 461, 109 
(1996) and references therein. 

[21] F. Berruto, G. Grignani and P. Sodano, Phys. Rev. D 62, 054510 (2000). 

[22] N. D. Mermin and H. Wagner, Phys. Rev. Lett. 17, 1133 (1966); S. Goleman, Gommun. 
Math. Phys. 31, 259 (1973). 

[23] A. R. Zhitnitsky, Phys. Lett. B 165, 405 (1985) [Sov. J. NucL Phys. 43, 999.1986 
YAFIA,43,1553 (1985)]. 

[24] E. Witten, NucL Phys. B 160, 57 (1979). 

[25] S. Goleman, Aspects of Symmetry, (Gambridge University Press, Gambridge, 1985). 
[26] P. H. Ginsparg and K. G. Wilson, Phys. Rev. D 25, 2649 (1982). 
[27] D. J. Gross and E. Witten, Phys. Rev. D 21, 446 (1980). 

[28] H. W. Hamber, E. Marinari, G. Parisi and G. Rebbi, Nucl. Phys. B 225, 475 (1983). 

[29] T. Eguchi and H. Kawai, Phys. Rev. Lett. 48, 1063 (1982). 

[30] H. Leutwyler and A. Smilga, Phys. Rev. D 46, 5607 (1992). 

[31] P. Hernandez, K. Jansen and L. Lellouch, Phys. Lett. B 469, 198 (1999); 

P. Hernandez, K. Jansen, L. Lellouch and H. Wittig, JHEP 0107, 018 (2001). 

[32] M. J. Teper, Phys. Rev. D 59, 014512 (1999); B. Lucini and M. Teper, JHEP 0106, 
050 (2001). 



17 



